Oh! 




Belle Prerpint 2006-10 
KEK Preprint 2006-4 



Measurements of branching fractions and q^ 
distributions for B — > niiy and B — > plu 
\b : Decays with B -^ D^*^iiy Decay Tagging 

o' 

(N; 

^' Belle Collaboration 

<: 

^ : T. HokuueS K. Abe^ K. Abe^™, I. Adachi^ H. Aihara^°, 

; Y. Asano^^ T. Aushev^, A. M. Bakich^J, V. Balagura^, 

> : E. BarberioP, M. Barbero'^, A. Bay™, I. Bedny^, K. Belous^ 

^ : U. Bitenc\ I. Bizjak\ S. Blyth^ A. Bondar^ A. Bozek^, 

§ : M. Bracko^'°'\ T. E. Browder^, P. Chang ^, A. Chen^ 

O : w. T. Chen^ Y. Choi^\ A. Chuvikov^'^, S. Cole^J, J. Dalseno?, 

O : M. Danilov^, M. Dash^S A. Dmtskoy^ S. Eidelman^ 

^ ! N. Gabyshev^, A. Garmash '^'^j T. Gershon^, G. Gokhroo^^, 

A. Gorisek\ H. Ha^, J. Haba^ K. Hara^ T. Hara^^ 
K. Hayasaka'i, H. Hayashii^, M. Hazumi'^, L. Hinz"^, 
Y. Hoshi^™, S. Hou^ W.-S. Hou'^, T. lijima^ K. Ikado<i, 
j>^ : A. Imoto^ K. Inami^i, A. Ishikawa''", R. Itoh"^, M. Iwasaki^°, 

^ '■ Y. Iwasaki^ H. Kakuno^°, J. H. Kang^^, P. Kapusta^, 

S. U. Kataoka'", H. Kawai^, T. Kawasaki "", H. R. Khan^P, 
H. J. Kim^, H. O. Kim^S K. Kinoshita^ Y. Kozakai^, 
P. Krizan'^'\ P. Krokovny^, R. Kulasiri*^, R. Kumar ^^, 
C. C. Kuo^ A. Kuzmin% Y.-J. Kwon^^, J. Lee^s, T. Lesiak^, 
J. Li^^ A. Limosani^ S.-W. Lin^, G. Majumder^^, F. Mandl^, 
T. Matsumoto^^ A. Matyja^, S. McOnie^J, W. Mitaroffs, 
H. Miyake^^, H. Miyata^, Y. Miyazaki'i, D. Mohapatra^*, 
I. Nakamura^, E. Nakano^, Z. Natkaniec^, S. Nishida*^, 
O. Nitoh^^ T. Nozaki^ S. Ogawa^^, T. Ohshima^, T. Okabe'i, 
S. OkunoJ, S. L. Olsen^, Y. Onuki^, P. Pakhlov^ 



Preprint submitted to Elsevier Science 7 February 2008 



ae 

5 



C. W. Park^\ H. Park^, L. S. Peak^J, R. PestotnikS 

L. E. Piilonen^*, Y. Sakai^ N. Sato^, N. Satoyama^^ 

K. Sayeed*^, T. Schietinger "^, O. Schneider"^, C. Schwanda^ 

A. J. Schwartz ^ K. Senyo^, M. E. SeviorP, M. Shapkin^ 

H. Shibuya^^, B. Shwartz^, A. Somov*^, R. Stamen '^j 

S. Stanic^, M. Staric\ H. Stoeck^J, A. Sugiyama^°, S. Suzuki 

S. Y. Suzuki^, O. Tajima'^, F. Takasaki^, K. Tamai*^, 

M. Tanaka°, G. N. Taylor?, Y. Teramoto^ X. C. Tian^^ 

T. Tsukamoto^ K. Ueno^, T. Uglov^, Y. Unno^ S. Uno^ 

P. UrquijoP, Y. Usov% G. Varner'^, K. E. VarvelPJ, S. Villa ^, 

C. C. Wang^, C. H. Wang*, M.-Z. Wang^, Y. Watanabe^P, 

E. Won'^, A. Yamaguchi ^^, Y. Yamashita'^, M. Yamauchi*^, 

J. Ying^^, L. M. Zhang ^*, Z. P. Zhang ^*, and D. Ziircher"^, 

^Budker Institute of Nuclear Physics, Novosibirsk, Russia 

Chiba University, Chiba, Japan 

^ University of Cincinnati, Cincinnati, OH, USA 

'^ University of Hawaii, Honolulu, HI, USA 

^High Energy Accelerator Research Organization (KEK), Tsukuba, Japan 

^Institute for High Energy Physics, Protvino, Russia 

^Institute of High Energy Physics, Vienna, Austria 

^Institute for Theoretical and Experimental Physics, Moscow, Russia 

V. Stefan Institute, Ljubljana, Slovenia 

^Kanagawa University, Yokohama, Japan 

^Korea University, Seoul, South Korea 

Kyungpook National University, Taegu, South Korea 

™ Swiss Federal Institute of Technology of Lausanne, EPFL, Lausanne, Switzerland 

" University of Ljubljana, Ljubljana, Slovenia 

° University of Maribor, Maribor, Slovenia 

P University of Melbourne, Victoria, Australia 

^Nagoya University, Nagoya, Japan 

^Nara Women's University, Nara, Japan 

^National Central University, Chung-li, Taiwan 

^National United University, Miao Li, Taiwan 

"Department of Physics, National Taiwan University, Taipei, Taiwan 

^H. Niewodniczanski Institute of Nuclear Physics, Krakow, Poland 

""Nippon Dental University, Niigata, Japan 



^Niigata University, Niigata, Japan 

^ Nova Gorica Polytechnic, Nova Gorica, Slovenia 

^ Osaka Gity University, Osaka, Japan 

'"^^ Osaka University, Osaka, Japan 

^^Panjab University, Chandigarh, India 

^Peking University, Beijing, PR China 

Princeton University, Princeton, NJ, USA 

^Saga University, Saga, Japan 

University of Science and Technology of China, Hefei, PR China 

^^ Seoul National University, Seoul, South Korea 

Shinshu University, Nagano, Japan 

^^ Sungkyunkwan University, Suwon, South Korea 

^^ University of Sydney, Sydney, NSW, Australia 

'^^ Tata Institute of Fundamental Research, Bombay, India 

^^Toho University, Funabashi, Japan 

^™ Tohoku Gakuin University, Tagajo, Japan 

'^^Tohoku University, Sendai, Japan 

^"Department of Physics, University of Tokyo, Tokyo, Japan 

^P Tokyo Institute of Technology, Tokyo, Japan 

'^Tokyo Metropolitan University, Tokyo, Japan 

^^' Tokyo University of Agriculture and Technology, Tokyo, Japan 

^ University of Tsukuba, Tsukuba, Japan 

Virginia Polytechnic Institute and State University, Blackshurg, VA, USA 

'^^ Yonsei University, Seoul, South Korea 



Abstract 

We report measurements of the charmless semileptonic decays B^ — > ■k~ / p~l~^i' 
and B^ —>■ -K^/p^i'^v, based on a sample of 2.75 x 10^ BB events collected at the 
T{4S) resonance with the Belle detector at the KEKB e~^e~ asymmetric collider. In 
this analysis, the accompanying B meson is reconstructed in the semileptonic mode 
B -^ D^*'iv, enabling detection of the signal modes with high purity. We measure 
the branching fractions B{B° -^ vr-^+i/) = (1.38 ± 0.19 ± 0.14 ± 0.03) x lO""^, 
B{B'^ -^ p~e+u) = (2.17 ± 0.54 ± 0.31 ± 0.08) x 10-^ 13{B+ -^ tt^I+u) = (0.77 ± 
0.14 ±0.08 ±0.00) X 10^^ and B{B+ -^ p°£+i/) = (1.33 ±0.23 ±0.17 ±0.05) x 10^^ 
where the errors are statistical, experimental systematic, and systematic due to 
form-factor uncertainties, respectively. For each mode we also present the partial 



branching fractions in three q^ intervals: (7^ < 8, 8 < (/^ < 16, and q^ > 16GeV^/c^. 
From our partial branching fractions for B — > Triv and recent results for the form 
factor from unquenched Lattice QCD calculations, we obtain values of the CKM 
matrix element iT^bl- 

Key words: Semileptonic, B decay, exclusive 
PACS: 12.15.Hh, 12.38.Gc, 13.25.Hw 



1 Introduction 



Exclusive B -^ X^iv decays proceed dominantly via a, h ^ uW~ tree pro- 
cess and can be used to determine | Vub \ , one of the smallest and least known 
elements of the Cabibbo-Kobayashi-Maskawa matrix [1]. However, the need 
to translate the observed rate to a \Vub\ value using model-dependent de- 
cay form-factors (FF) has resulted in large theoretical uncertainties. The re- 
cent release of FF results for B -^ -irii/ calculated by unquenched Lattice 
QCD (LQCD) [2,3] makes possible the first model-independent determina- 
tion of \Vub\- Since LQCD results are available only in the high q^ region 
(> 16 GeV^/c^), a clean measurement of the partial B -^ irlu branching frac- 
tion in the same high g^ region is needed. 

There have been several measurements in the past by CLEO, BaBar and Belle 
for the B — > niu, piu, rjlv and ujIi' modes [4,5,6,7,8,9]. The analyses in these 
measurements utilize the method, originally developed by CLEO, where the 
B decays are reconstructed by inferring the undetected neutrino mass from 
missing energy and momentum ( 'V- reconstruction method") [4]. In the B- 
factory era, we will improve the statistical precision by simply applying the 
i^-reconstruction method using a large amount of data. However, the poor 
signal-to- noise ratio will limit the systematic uncertainty of the measurement. 

In this paper we present measurements of B^ -^ tx~ j p~t^v and B^ — * 
'K^ I p^t^v decays using B — ^ D'^*^lv decay tagging. We reconstruct the en- 
tire decay chain from the T(4S'), T(4S') -^ -Bsig-Btag, -Bsig -^ n/ pdu and 
-Stag -^ D^*^!/ with several D^*^ sub-modes. The back-to-back correlation 
of the two B mesons in the T(45') rest frame allows us to constrain the kine- 
matics of the double semileptonic decay. The signal is reconstructed in four 
modes, B^ -^ 7r~ /p~i~^i' and B~^ —>■ vr^/p^f+z/. Yields and branching frac- 
tions are extracted from a simultaneous fit of the B^ and B~^ samples in three 
intervals of g^, accounting for cross- feed between modes as well as other back- 
grounds. We have applied this method to B ^ ir/piu decays for the first time, 
and have succeeded in reconstructing these decays with significantly improved 
signal-to-noise ratios compared to the //-reconstruction method. Inclusion of 
charge conjugate decays is implied throughout this paper. 



2 Data Set and Experiment 



The analysis is based on data recorded with the Belle detector at the KEKB 
collider operating at the center-of-mass (cm.) energy for the T{4S) reso- 
nance [10]. The T(4S') dataset that is used corresponds to an integrated lu- 
minosity of 253 fb~^ and contains 2.75 x 10® BB events. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of 
a silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an ar- 
ray of aerogel threshold Cerenkov counters (ACC), a barrel-like arrangement of 
time-of-flight scintillation counters (TOF), and an electromagnetic calorime- 
ter comprised of CsI(Tl) crystals (ECL) located inside a super-conducting 
solenoid coil that provides a 1.5 T magnetic field. An iron fiux-return located 
outside of the coil is instrumented to detect K^ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [11]. Two inner detector 
configurations were used. A 2.0 cm beam pipe and a 3-layer silicon vertex 
detector was used for the first sample of 152 x 10^ BB pairs, while a 1.5 cm 
beam pipe, a 4-layer silicon detector, and a small-cell inner drift chamber were 
used to record the remaining 123 million BB pairs [12]. 

A detailed Monte Carlo (MC) simulation, which fully describes the detector 
geometry and response and is based on GEANT [13], is applied to estimate 
the signal detection efficiency and to study the background. To examine the 
FF dependence, MC samples for the B -^ Ti/piu signal decays are generated 
with different form-factor models: a quark model (ISGW II [14]), light cone 
sum rules (LCSR for B — > -Kiu [15] and B — >■ piu [16]) and quenched lattice 
QCD (UKQCD [17]). We also use unquenched lattice QCD (FNAL [2] and 
HPQCD [3]) for B -^ Triu and a relativistic quark model (Melikhov [18]) 
for B -^ piu. To model the cross-feed from other B -^ X^iu decays, MC 
samples are generated with the ISGW II model for the resonant components 
{B —>■ 7rii> and B —>■ piu components are excluded in this sample) and the 
DeFazio-Neubert model [19] for the non-resonant component. To model the 
BB and continuum backgrounds, large generic BB and qq Monte Carlo (based 
on Evtgen [20]) samples are used. 



3 Event Reconstruction and Selection 



Charged particle tracks are reconstructed from hits in the SVD and the CDC. 
They are required to satisfy track quality cuts based on their impact parame- 
ters relative to the measured interaction point (IP) of the two beams. Charged 
kaons are identified by combining information on ionization loss (dE/dx) in the 
CDC, Cherenkov light yields in the ACC and time-of-flight measured by the 



TOF system. For the nominal requirement, the kaon identification efficiency 
is approximately 88% and the rate for misidentification of pions as kaons is 
about 8%. Hadron tracks that are not identified as kaons are treated as pi- 
ons. Tracks satisfying the lepton identification criteria, as described below, are 
removed from consideration. 

Neutral pions are reconstructed using 7 pairs with an invariant mass between 
117 and 150 MeV/c^. Each 7 is required to have a minimum energy of 50 MeV. 
Kg mesons are reconstructed using pairs of tracks that are consistent with 
having a common vertex and that have an invariant mass within ±12 MeV/c^ 
of the known Kg mass. 

Electron identification is based on a combination of dE/dx in the CDC, the 
response of the ACC, shower shape in the ECL and the ratio of energy de- 
posit in the ECL to the momentum measured by the tracking system. Muons 
are identified by their signals in the KLM resistive plate counters, which are 
interleaved with the iron of the solenoid return yoke. The lepton identification 
efficiencies are estimated to be about 90% for both electrons and muons in 
the momentum region above 1.2 GeV/c, where leptons from prompt B decays 
dominate. The hadron misidentification rate is measured using reconstructed 
K^ -^ TT^TT^ and found to be less than 0.2% for electrons and 1.5% for muons 
in the same momentum region. 

For the reconstruction of -Btag -^ D^*^^!?, the lepton candidate is required to 
have the correct sign charge with respect to the D meson flavor and a labora- 
tory momentum (p^"*) greater than 1.0 GeV/c. The D meson candidates are re- 
constructed by using seven decay modes of D~^: D'^ -^ K^n^ii^, K~'7t'^7i~^7t^, 
Kg7r+, K^TT+vr", /f^TT+TT+TT-, K+K'^g, K+K-TT+; and ten decay modes of D°: 
D^ -^ K-TT+, A'-7r+7r°, K'tt+tt+tt-, K^^ti^ , i^^TT+vr-, K^ti+ti-ti^, K- ti+ ti+ ti' -k^ 
K~^K~ , KgK~^K~ , KgK~7i~^. The candidates are required to have an invari- 
ant mass m£) within ±2.5cr [a is a standard deviation) of the nominal D mass, 
where the mass resolution a is dependent on the decay mode. D* mesons are 
reconstructed in the modes D*~^ — > D^tt"*", D~^7r^ and D*'^ — *> D^ir^ by combin- 
ing a D meson candidate and a charged or neutral pion. Each D* candidate 
is required to have a mass difference Am = mj^^ — rrijj within ±2. 5a of the 
nominal values. 

For the reconstruction of -Bgig -^ Xud-V-, the lepton candidate is required to 
have the right sign charge with respect to the X^ system and p^"* greater than 
0.8 GeV/c. The X^ system may consist of one pion or two pions ( Nt^+ = 1 or 
Nt,+ = Nt,o = 1 for a i?o tag and A^'^o = 1 or Nt,+ = N^,- = 1 for a B~ tag). 
The event is required to have no additional charged tracks or 7r° candidates. 
We also require that the residual energy from neutral clusters be less than 
0.15 GeV (-Eneut < 0.15 GeV). The two leptons on the tag and the signal sides 
are required to have opposite charge. The loss of signal due to 5° — 5° mixing 



is estimated by MC simulation. 

We then impose a constraint based on the kinematics of the double semilep- 
tonic decay in the T(4S') rest frame. In the semileptonic decay on each side, 
Bi(2) -^ yi{2)^ iXi = D^*H and Y2 = X^i), the angle between the Si(2) 
meson and the detected ^1(2) system Obj^,^) ^^ calculated from the relation, 



p* 



{Pe 



P^f 



(P*: 4-momentum vector) and the known p*^ (the 



absolute momentum of the mother B meson). This means that the -81(2) di- 
rection is constrained on the surface of a cone defined with the anerle 9%, 

^ -°1(2) 

around the direction of the 11(2) system, as shown graphically in Fig. 1. The 
back-to-back relation of the two B meson directions then implies that the 
real B direction is on the intersection of the two cones when one of the B 
systems is spatially inverted. Denoting 912 the angle between the py^ and 
— Py25 the B directional vector ub = {xb.Ub^zb) is given by, zb = cos6'|j^, 
i/B = {cos9*B^ - cos9*B^cos9l2)/sm9l2, and 



Xb 



— 97;— (cos^6'r, 
sin^^i^ ' 
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B2 



2cos6'^j cos^Ij^ cos^^2 ) 



:r 



'12 



with the coordinate definition in Fig. 1, where the pyi and pyg ^^^ aligned 
along the 2;-axis and in the y — z plane, respectively. If the hypothesis of the 
double semileptonic decay is correct and all the decay products are detected 
except for the two neutrinos, x^ must range from to 1. Events passing a 
rather loose cut x^ > —2.0 are used for signal extraction at a later stage of 
the analysis. 
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Fig. 1. Kinematics of the double semileptonic decay. 
Since the direction of the B meson is not uniquely determined, we calculate 

V beam 



2 2 

q as q 



E*xy -Pk^ using the beam energy (^beam), energy 



(Ex^) and momentum {p*x^) of the Xu system and neglecting the momentum 
of the B meson in the cm. system. The signal Monte Carlo simulation finds 
that the g^ resolution depends on the reconstructed g^ and is in the range 
0.32-0.95 GeV7c2. 

According to Monte Carlo simulation, the largest backgrounds originate from 
B —>■ Xciiy and non-signal B —>■ X^^v decays, where some particles escape 
detection. There are sizable contributions from cross talk between the B^ and 
B'^ tags. The contribution from qq processes is found to be negligible. 

For events selected as described above, the signal MC simulation indicates 
that the total detection efficiency (etotai), averaged over electron and muon 
channels, is 1.98 x lO'^ for 5° -^ ir-i+u and 0.76 x lO'^ for S° -^ p'i^i^, 
1.49 X 10-3 for B+ -^ n^i+u and 1.78 x lO'^ for B+ -^ p^i+u assuming 
the LCSR FF model. Here, etotai is defined with respect to the number of 
BB pairs, where one B decays into the signal mode, and includes the loss of 
signal due to B^ — B^ mixing. Because of the loose lepton momentum cut (p^"^ 
> 0.8 GeV/c), the variation of efficiency with different FF models is relatively 
small. Table 1 gives the matrix e{qrec.yQtTne)y the efficiency for a signal event 
generated with true g^ in the bin g^^.^^ to be reconstructed in the bin g^^^ . 

Table 1 
Detection efficiency matrix e(9rcc. ' 9?rue) based on the LCSR model in units of 10"^. 





4jGeVVc2) 


2 
Qtrue 


TT-i+U 


p~l+u 


vrO^+z. 


p^e+i^ 


(GeVVc2) 


< 8 8 - 16 > 16 


< 8 8 - 16 > 16 


< 8 8- 16 > 16 


< 8 8-16 > 16 


<8 
8-16 
> 16 


1.71 0.21 0.00 
0.05 1.82 0.24 
0.00 0.03 1.89 


0.59 0.03 0.00 
0.07 0.65 0.05 
0.02 0.13 0.81 


1.27 0.07 0.00 
0.10 1.43 0.06 
0.01 0.09 1.45 


1.50 0.08 0.00 
0.10 1.71 0.13 
0.01 0.08 1.82 



To check the validity of the method, we apply the procedure described above 



to reconstruct 5°g -^ D*-i+u followed by D*' -^ D^n^ , D^ -^ K+tx- for 
a 5° tag and S+g -^ D*^i+u followed by D*^ -* D^tx^, D^ -^ K+ti^ for a 
B~ tag, with the same requirement on the tagging side. Figures 2-a) and c) 
show the Mxtttt distributions that are obtained in data and expected from 
MC. As a resuh, we obtained 224.7 ± 15.4 (295.9 ± 17.6) B^ {B~) tagged 
events. These values are in good agreement with expected values 224.5 ± 9.5 
(288.6 ± 11.7) calculated from the branching fractions B{B^ —>■ D*~i~^i>), 
B{D*-^^^ -^ DOvr-W) and B{D^ -^ i^+vr") in [21] and efficiencies obtained 
from MC. Here, we use B{B+ -^ D*^(i+u) calculated from B{B^ -^ D*-(i+u) 



and the hftime ratio [21]; B{B^ -^ D 



*0/ 



V) 



B^B'^ -^ D*~(i+u) X {tb+/tbo] 



The ratio of the reconstructed to expected value, R = 1.00±0.08±0.05 (1.03± 
0.07 ± 0.05) where the first error is statistical error and the second is due to 
the uncertainty of the branching fractions from [21], is consistent with unity. 
Figures 2-b) and d) show a comparison of the reconstructed x% distribution 



in the above data samples with MC simulation. Data and MC are in good 
agreement. 
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Fig. 2. Reconstructed M^Kmr) distribution (a) and x^ distribution (b) for the 
B^ -^ D*~i^u calibration decay, (c) and (d) are for the B~^ —>■ D*^!t'^v decay; 
points with error bars are data and the histogram is the signal MC. 



4 Extraction of Branching Fractions 



The B^ — > TT~ / p~l'^v and B^ — > -n^ j p^t^v signals are extracted using binned 
maximum likelihood fits to the two-dimensional {x\^ Mx) distribution, where 
Mx is the nominal pion mass for B -^ ni'^v candidates and the invariant 
mass of two pions for B —^ pl^v candidates. The fit includes seven compo- 
nents: the four signal modes and the other B^ -^ X~i^v and B^ -^ Xl^i'^u 
backgrounds, the background from BB events containing no -B -^ X„£z/. The 
PDF (probability density function) for each fit component is determined from 
MC simulation. The vr/p signal events exhibit characteristic behavior in both 
the x% and Mx distributions; other B —>■ Xui'^u events exhibit a weak peak- 



ing structure in x^ but a broad distribution in Mx'-, the BB background has 
a relatively fiat distribution in x\ and a broad structure in Mx- The PDFs 
in (x^,Mx) for each of the seven fit components are obtained from MC for 
both B^ and B~ tag candidates. We then fit the two {x\,Mx) distributions 
for both B^ and B~ tags simultaneously; The fitting is constrained so that the 
sum of the deduced branching fractions for B -^ ni'^h', B -^ pi~^u and B -^ 
other Xui'^i' is equal to the total inclusive branching fraction B{B -^ Xulv) = 
(0.25±0.06)% [22]. Figure 3 shows the projections on Mx and x"^ of the fitting 
result for data in the entire g^ region. The extracted yields for the signal com- 
ponents are N{B^ -^ n^i+u) = 155.8 ± 20.0, N{B^ - 
N{B+ -^ TT^i+u) = 69.0 ± 11.4 and N{B+ -^ p^i+u) 



p-e+u) =92.9 ±19.4, 
135.4 ±24.8, with the 



LCSR model used for the four signal PDFs. 
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Fig. 3. Projected Mx distribution for x^ > —2 (left) and x^ distributions for 
the mass region of vr {Mx < O.lSGeV/c^, middle) and p (0.6 < Mx < l.OGeV/c^, 
right) in all q^ region; points are data. Histogram components are B^ -^ 'K~i^u (red 
narrow 135° hatching), B^ — > p~l'^v (red wide 45° hatching), other X^t^v from 
B^ (red cross-hatching) and B^ —> n^i'^iy (blue narrow 45° hatching), B'^ -^ p^i'^v 
(blue wide 135° hatching), other Xui'^v from B^ (blue cross-hatching) and BB 
background (green horizontal hatching). 

Figure 4 shows projections of the data, separated into three q^ bins, g^ < 8 
GeV^/c^, 8 < g^ < 16 GeV^/c^ and q^ > 16 GeV^/c^. Here the normalizations 
of the other B -^ X^^v and the BB background components are fixed to 
those obtained in the above fitting for the entire q^ region. The extracted 



10 



numbers of events for low/ medium/ high q^ bins are, N{B^ -^ ir i^u) = 
64.8 ± 11.9 / 63.2 ± 12.4 / 40.6 ± 11.3, N{B^ -^ p-i+u) = 22.1 ± 8.0 / 53.2 ± 
13.5 / 30.9 ± 16.0, N{B+ -^ ttH+u) = 18.1 ± 5.1 / 34.5 ± 8.3 / 18.6 ± 6.5 
and N{B+ -^ p°£+z/) = 47.2 ±11.2 / 68.3 ± 16.5 / 32.5 ± 12.3. Table 2 
summarizes the extracted branching fractions. The branching fractions are 
calculated for each signal FF-model, where we take the average for cross-feed 
FF-models. The results are unfolded using the efficiency matrix e(g^gj. , g^uc) 
for the three g^ intervals prepared for each signal FF-model. We calculate the 
total branching fraction by taking the sum of the partial branching fractions 
in the three g^ intervals. 
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Fig. 4. Projected Mx distribution for x^ > —2 in each q^ region; points are 
data. Histogram components are B^ -^ ■k~(.'^v (red narrow 135° hatching), 
B^ —f p~l'^v (red wide 45° hatching), other Xu^'^v from B^ (red cross-hatching) and 
B^ -^ -n^l'^iy (blue narrow 45° hatching), B'^ —> p^i'^u (blue wide 135° hatching), 
other Xui'^i' from B~^ (blue cross-hatching) and BB background (green horizontal 
hatching) . 



5 Systematic Errors 



Tables 3 and 4 summarize the experimental systematic errors on the branching 
fractions. The experimental systematic errors can be categorized as originat- 
ing from uncertainties in the signal reconstruction efficiency, the background 
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Table 2 

Extracted branching fractions for each signal mode with different FF models in units 
of 10~^: the total branching fraction and the partial branching fractions in three q^ 
intervals, x^ /n and the associated probability for this x^ indicate the quality of the 
fit for the FF shape to the observed q^ distribution. 



Mode 


Model 


^total 


B<8 


^8-16 


13>w 


X^/n 


Prob. 


■K-(.+v 


LCSR 


1.40 ±0.19 


0.52 ±0.11 


0.51 ±0.11 


0.36 ±0.10 


0.4/2 


0.81 




ISGW II 


1.36 ±0.19 


0.52 ±0.11 


0.48 ±0.11 


0.36 ±0.10 


3.6/2 


0.17 




UKQCD 


1.39 ±0.19 


0.53 ±0.11 


0.50 ±0.11 


0.36 ±0.10 


0.2/2 


0.89 




FNAL 


1.39 ±0.19 


0.52 ±0.11 


0.51 ±0.11 


0.36 ±0.10 


0.3/2 


0.86 




HPQCD 


1.39 ±0.19 


0.52 ±0.11 


0.50 ±0.11 


0.36 ±0.10 


0.5/2 


0.79 




Average 


1.38 ±0.19 


0.52 ±0.11 


0.50 ±0.11 


0.36 ±0.10 


- 


- 


p-i+u 


LCSR 


2.21 ±0.54 


0.48 ±0.19 


1.18 ±0.33 


0.55 ±0.39 


2.4/2 


0.43 




ISGW II 


2.09 ±0.54 


0.44 ±0.17 


1.14 ±0.32 


0.51 ±0.39 


0.5/2 


0.78 




UKQCD 


2.19 ±0.54 


0.48 ±0.19 


1.20 ±0.33 


0.51 ±0.38 


1.7/2 


0.43 




Melikhov 


2.20 ±0.54 


0.48 ±0.19 


1.22 ±0.34 


0.50 ±0.37 


1.8/2 


0.41 




Average 


2.17 ±0.54 


0.47 ±0.19 


1.19 ±0.33 


0.52 ±0.38 


- 


- 


71^1+ U 


LCSR 


0.77 ±0.14 


0.19 ±0.06 


0.39 ±0.10 


0.20 ±0.08 


2.9/2 


0.24 




ISGW II 


0.77 ±0.14 


0.19 ±0.06 


0.39 ±0.10 


0.20 ±0.08 


7.8/2 


0.02 




UKQCD 


0.77 ±0.14 


0.19 ±0.06 


0.38 ±0.10 


0.20 ±0.08 


2.5/2 


0.28 




FNAL 


0.77 ±0.14 


0.19 ±0.06 


0.39 ±0.10 


0.20 ±0.08 


2.8/2 


0.25 




HPQCD 


0.77 ±0.14 


0.19 ±0.06 


0.39 ±0.10 


0.20 ±0.08 


4.8/2 


0.09 




Average 


0.77 ±0.14 


0.19 ±0.06 


0.39 ±0.10 


0.20 ±0.08 


- 


- 


p°i+l^ 


LCSR 


1.36 ±0.23 


0.43 ±0.11 


0.63 ±0.16 


0.29 ±0.12 


0.9/2 


0.64 




ISGW II 


1.29 ±0.22 


0.45 ±0.11 


0.57 ±0.15 


0.27 ±0.12 


1.4/2 


0.50 




UKQCD 


1.34 ±0.23 


0.45 ±0.11 


0.62 ±0.16 


0.28 ±0.12 


0.2/2 


0.91 




Melikhov 


1.35 ±0.23 


0.44 ±0.11 


0.63 ±0.17 


0.28 ±0.12 


0.2/2 


0.92 




Average 


1.33 ±0.23 


0.44 ±0.11 


0.61 ±0.16 


0.28 ±0.12 


- 


- 



estimation, and the normalization. The total experimental systematic error is 
the quadratic sum of all individual ones. We also consider the systematic error 
due to the dependence on the FF model. 

The effect from the uncertainty on the signal reconstruction efficiency is eval- 
uated based on the efficiency calibration with the B, 



Slg 



D*(.'^u sample, 

discussed above. The error is taken to be that on the ratio of observed to ex- 
pected number of the calibration signals (9.3% for B^ -^ tx~ /p~t^v, 9.2% for 
B — »• t:^ / p^t^v). This gives the largest contribution to the systematic error. 
Note that this error is dominated by the statistics of the calibration signals, as 
explained above. Therefore, accumulation of additional integrated luminosity 
in the future will help to reduce this uncertainty. We also include residual 
errors for the reconstruction of the signal side: 1% and 2% for the detection 
of each charged and neutral pion, respectively, and 2% for the charged pion 
selection and 2.1% for the lepton selection. 
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The systematic error due to the uncertainty on the inclusive branching frac- 
tion B{B -^ X„£z/), which is used to constrain B — *> Xui'^i' background, is 
estimated by varying this parameter by its ±lo" error. The uncertainty in 
the BB background shape after our pion multiphcity selection requirements 
[N^+ = 1 or N^+ = N^o = 1 for a B^ tag and A^^o = 1 or N^+ = N^- = 1 for 
a B' tag) is studied in the simulation by randomly removing charged tracks 
and 7r° according to the error in detection efficiency (1% for a charged track, 
2% for 7r°), and also by reassigning identified charged kaons as pions accord- 
ing to the uncertainty in the kaon identification efficiency (2%). The resultant 
changes in the extracted branching fractions are assigned as systematic errors. 
We find a significant uncertainty in the high g^ region (g^ > 16 GeV^/c^) for 
B — >■ pi'^u due to the poor signal-to-noise ratio. We also vary the fraction of 
B — s> D**lv decays in the BB background MC by the error quoted in [21] to 
test the B — > Xcif model dependence in the BB background shape. To assess 
the uncertainty due to the production rate of K^, we vary the production rate 
in the MC simulation by the uncertainty in the inclusive branching fraction 
for B ^ K° X quoted in [21]. 

For the normalization, we consider the uncertainty in the number of B^B^ and 
B+B- pairs: the ratio of B+B' to B^B^ pairs (U/fo), U/fo = 1.029 ±0.035 
[23] , the mixing parameter (xd), Xd = 0.186 ± 0.004 [21], and the measured 
number of BB pairs {N^j^, 1.1%). 

Table 3 

Summary of systematic errors (%) for B{B^ —> i:^ j p^ t^v\ 







B^^ 


vr-^+i/ 




B^- 


> p^t 


v 




l' 


interval (GeVVc^) 


5^ 


interval (GeVVc^) 


Source 


g2 < 8 8 - 16 > 16 


< 16 


all 


g2 <8 


8 - 16 > 16 


< 16 


all 


Tracking efficiency 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


TT^ reconstruction 


- 


- 


- 


- 


- 


2 


2 


2 


2 


2 


Lepton identification 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


Kaon identification 


2 


2 


2 


2 


2 


2 


2 


2 


2 


2 


D*iv calibration 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


9.3 


BriXutv) in the fitting 


0.8 


2.4 


1.8 


1.6 


1.4 


4.8 


3.9 


23.8 


1.9 


7.1 


BB background shape 


1.1 


2.2 


2.8 


1.2 


1.3 


3.8 


2.9 


17.0 


3.0 


6.1 


Br{D**£u) 


1.0 


1.5 


0.2 


1.2 


0.9 


0.5 


0.3 


2.5 


0.3 


0.8 


K^ production rate 


0.1 


0.3 


0.4 


0.2 


0.3 


1.0 


0.7 


2.9 


0.8 


1.3 


^BB 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


f+/fo 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


Xd 


2.2 


2.2 


2.2 


2.2 


2.2 


2.2 


2.2 


2.2 


2.2 


2.2 


exp. total 


10.4 


10.9 


10.8 


10.5 


10.5 


12.1 


11.5 


31.3 


11.1 


14.1 


FF for signal 


0.7 


3.8 


0.9 


2.2 


1.8 


6.1 


3.5 


6.8 


4.3 


3.6 


FF for cross-feed 


1.8 


2.1 


1.5 


1.9 


1.4 


0.5 


0.7 


2.4 


0.6 


1.0 


FF total 


1.9 


4.3 


1.7 


2.9 


2.3 


6.1 


3.6 


7.2 


4.3 


3.7 
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Table 4 

Summary of systematic errors (%) for B{B'^ -^ -n^ / p^ (.'^ v) . 







B+ - 


> TT^i+V 




B+ - 


^pO^-t 


v 




q' 


interva 


I (GeVVc2) 


2 

q 


interva 


( GeVVc2) 


Source 


q^<8 


8-16 > 16 


< 16 


all 


q^<8 


8 - 16 > 16 


< 16 


all 


Tracking efficiency 


- 


- 


- 


- 


- 


2 


2 


2 


2 


2 


■K^ reconstruction 


2 


2 


2 


2 


2 


- 


- 


- 


- 


- 


Lepton identification 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


2.1 


Kaon identification 


- 


- 


- 


- 


- 


4 


4 


4 


4 


4 


D*lu calibration 


9.2 


9.2 


9.2 


9.2 


9.2 


9.2 


9.2 


9.2 


9.2 


9.2 


Br{XJv) in the fitting 


0.2 


3.1 


3.0 


2.1 


1.2 


2.0 


3.7 


20.0 


3.0 


6.6 


BB background shape 


1.9 


5.5 


2.7 


4.3 


3.7 


5.3 


4.3 


16.3 


1.5 


2.8 


Br{D**£u) 


1.3 


0.8 


0.8 


0.9 


0.9 


0.2 


1.6 


3.0 


0.9 


1.4 


K^ production rate 


0.3 


1.1 


0.6 


0.8 


0.8 


0.3 


0.2 


1.9 


0.1 


0.4 


^BB 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


1.1 


f+/fo 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


1.7 


exp. total 


10.1 


11.8 


10.7 


11.0 


10.7 


12.1 


12.2 


28.1 


11.2 


12.9 


FF for signal 


1.2 


0.5 


1.3 


0.3 


0.2 


2.1 


7.1 


3.9 


3.7 


3.5 


FF for cross-feed 


0.7 


0.8 


0.6 


0.8 


0.6 


3.3 


1.1 


1.0 


1.5 


1.2 


FF total 


1.4 


0.9 


1.4 


0.8 


0.6 


3.9 


7.2 


4.0 


4.0 


3.7 



The dependence of the extracted branching fractions on the FF model has 
been studied by repeating the above fitting procedure with various FF models 
for the signal mode and also for the cross-feed mode {B -^ iriu ^^ B ^ 
piu). We consider the models listed in Table 2. For the extracted B{B — >■ 
TT'f+z/ (tt^^+z/)), the standard deviation among the models is < 1.7 (0.9)% 
for B -^ TTi+u and < 1.9 (0.5)% for B -^ pi+u. For 13{B -^ p^i+u (p°£+z/)), 
the standard deviation is < 2.9 (3.6)% for B -^ pi+u and < 1.0 (1.3)% for 
B -^ ni^u. The total error due to FF model dependence is the quadratic sum 
of the maximum variations with the signal and cross-feed FF models. 



6 Results 



Table 5 summarizes our measurements of the total and partial branching frac- 
tions for the four signal modes. Each branching fraction is obtained by tak- 
ing the simple average of the values obtained from the FF models shown in 
Table 2. The errors shown in the table are statistical, experimental system- 
atic, and model dependence due to form-factor uncertainties. The obtained 
branching fractions for B^ -^ 7r~/p"£+i/ are consistent with the existing mea- 
surements by CLEO [6] and BaBar [9]. The overall uncertainty on our result 
for B^ -^ 7r~ i'^i/ (17%) is comparable to those on CLEO and BaBar results 
based on //-reconstruction. Our results for B^ -^ p~i,^v have the smallest 
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uncertainty. 

Table 5 

Summary of the obtained branching fractions. The errors are statistical, experimen- 
tal systematic, and systematic due to form-factor uncertainties. 



Modes q^ region (GeV^/c^) 


Branching fraction (x 10 ^) 


B^ -^ TT-i+u Total 

> 16 
< 16 


1.38 ±0.19 ±0.14 ±0.03 
0.36 ±0.10 ±0.04 ±0.01 
1.02 ±0.16 ±0.11 ±0.03 


B+ -^ TT^i+i, Total 

> 16 
< 16 


0.77 ±0.14 ±0.08 ±0.00 
0.20 ±0.08 ±0.02 ±0.00 
0.57 ±0.12 ±0.06 ±0.00 


B° ^ p-l+v Total 


2.17 ±0.54 ±0.31 ±0.08 


B+ -^ p^i+u Total 


1.33 ±0.23 ±0.17 ±0.05 



Figure 5 presents the measured q^ distributions for each signal mode, overlaid 
with the best fits of FF shapes to the data. To be self-consistent, the shape of 
a particular FF model is fit to the q^ distribution extracted with the same FF 
model. The quality of the fit in terms of x^ ^ind the probability of x^? shown 
in Table 2, may provide one way to discriminate among the models. From our 
results, the ISGW II model is disfavored for B — ^ nl^v. 

In this work, the 5° -^ Tr~i+u/B+ -^ ir^i+u and fi° -^ p-l+v/B+ -^ p^i+u 
signals are extracted separately, which allows us to test the isospin relations. 
From the obtained branching fractions and the B meson lifetimes in [21], the 
ratios of decay rates are found to be. 



r(E° ^ p~i+u) 

V{B+ -^ pOi+u) 



;i.92± 0.43 ±0.28), 
(1.74 ±0.53 ±0.33), 



(3) 



where the first and second errors are statistical and systematic errors, re- 
spectively. Both ratios are found to be consistent with the isospin relations; 
1(5° ^ Tr~{p~)i+iy) = 2T{B+ -^ 7rO(pO)£+i/). 

The obtained branching fractions in Table 5 can be used to extract IV^^I using 
the relation, 



IK 



ub\ 



B{B -^ TTJ+l^) 



(4) 



where Tthy is the form-factor normalization, predicted from theories. In Ta- 
ble 5, we list the partial branching fractions for B -^ iri'^i' decays in the q^ 
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Fig. 5. Extracted q^ distribution. Data points are shown for different FF models 
used to estimate the detection efficiency. Lines are for the best fit of the FF shapes 
to the obtained q^ distribution. 

region above 16 GeV^/c^, where the LQCD calculations are most reliable. The 
table provides also the results in the region below 16 GeV^/c^, so that one can 
deduce \Vub\ based on other approaches such as LCSR calculations [15,16]. 

In this paper we calculate \Vub\ based on the B —>■ ■nl'^v data in the high q^ 
region and the form factor predicted by recent unquenched LQCD calcula- 
tions. Their predictions {Tthy) for the q^ > 16GeV^/c^ region are Vthy{B'^ -^ 
Ti-i+u) = 1.83 ± 0.50 ps^i (FNAL) [2] and tthy{B° -^ 7r~i+u) = 1.46 ± 0.35 
ps-i (HPQCD) [3]. We use tbo = 1.532 ±0.009 ps and tb+ = 1.638 ±0.011 ps 
[21], and we use isospin symmetry to relate Tt^y for B^ -^ vr" and B^ — *> 7r° 
transitions. The results for B^ -^ tt'I^u and B^ — > tt'^^+z/ are then averaged, 
weighted by their respective statistical errors. 
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Table 6 

Summary of \Vub\ obtained from the B -^ Tri^i^ data in the q^ > 16GeV^/c^ 
region. The first and second errors are experimental statistical and systematic errors, 
respectively. The third error stems from the error on T^f^y quoted by the LQCD 
authors, 



Theory 


^thyips 1) 


Mode 


|K6|(xlO-3) 


FNAL 


1.83 ±0.50 




3.59 ± 0.51 ±0.20to:4i 
3.63 ± 0.70 ± 0.20to:4i 

3.60 ± 0.41 ±0.20to:4i 


HPQCD 


1.46 ± 0.35 




4.02 ± 0.57 ±0.22to:4? 
4.06 ± 0.78 ±0.22to:4? 

4.03 ± 0.46 ± 0.22to:4? 



Table 6 summarizes the results, where the first and second errors are the 
experimental statistical and systematic errors, respectively. The third error is 
based on the error on Tthy quoted by the LQCD authors. These theoretical 
errors are asymmetric because we assign them by taking the variation in \Vub\ 
when Tthy is varied by the quoted errors. The values are in agreement with 
those from inclusive B -^ X^^u decays [24]. 

To summarize, we have measured the branching fractions of the decays B — > 
TT-^i/ and B —>■ piu in 2.75 x 10® BB events using a method which tags one B in 
the mode B -^ D^*^(.v. Our results are consistent with previous measurements, 
and their precision is comparable to that of results from other experiments. 
The ratios of results for neutral and charged B meson modes are found to be 
consistent with isospin. The partial rates are measured in three bins of q^ and 
compared with distributions predicted by several theories. From the rate in 
the region q^ > 16 GeV^/c^ and recent results from LQCD calculations, we 
extract \Vub\'- 



,n-i+u+7rOe+u _ /Q en 4- n /11 -u n on+o.62^ 



l^-^'lKie GevVc^) = (3-60 ± 0.41 ± 0.20+^1) x 10-'(FNAL LQCD), (5) 

I(g2>l6 GcV2/c2) 



iV-blJ^mclvVc^) = (4-03 ± 0.46 ± 0.22+°^?) x 10-=^(HPQCD LQCD).(6) 



The experimental precision on these values is 13%, currently dominated by the 
statistical error of 11%. By accumulating more integrated luminosity, a mea- 
surement with errors below 10% is feasible. With improvements to unquenched 
LQCD calculations, the present method may provide a precise determination 

of \Vub\- 
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